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ABSTRACT 

Computational  analysis  and  two  dimensional  ten¬ 
sile  tests  on  single  motor  grain  fins  suggest  that 
cracks  in  fin  tips  are  most  likely  to  originate  at  the  co- 
'  alescence  of  the  fin  end  tip  radius  with  a  small  radius 
emanating  from  the  side  of  the  fin.  Prior  studies  have 
indicated  that  under  internal  pressure,  cracks  on  the 
fin  axis  are  subject  to  similar  stress  peaks  and  may 
grow  more  readily  than  the  former  types  due  to  aii 
absence  of  shear  modes.  The  present  study  focuses 
upon  two  types  of  cracks  emanating  from  the  former 
location  called  “off-axis”  cracks  and  attempts  to  dif¬ 
ferentiate  from  the  two  types  by  their  paths  and  SIF 
values,  determined  by  the  frozen  stress  photoelastic 
method. 

INTRODUCTION 

Tension  tests  on  pie  shaped  sections  of  a  generic 
motor  grain  model  containing  a  single  fin  at  AFRL 
revealed  that  cracks  initiated  at  the  confluence  of 
the  central  fin  tip  radius  and  the  edge  radius.  A  pre¬ 
liminary  photoelastic  study  of  an  uncracked  fin  has 
shown  (Smith,  Constantinescu  and  Liu,  2001)  that 
the  peak  stress  on  the  fin  axis  tip  is  nearly  as  great  as 
at  the  off-axis  position.  In  order  to  study  the  growth 
of  cracks  from  the  two  above  noted  positions  and 
measure  the  stress  intensity  factor  (SIF’s)  at  various 
locations  along  the  crack  borders,  a  series  of  frozen 
stress  (Appendix  A)  photoelastic  experiments  were 
carried  out  on  generic  six  finned  motor  grain  models 
containing  two  cracks  in  each  model  (Fig.  1).  The 


cracks  were  assumed  to  grow  independently  of  each 
other. 

THE  EXPERIMENTS 
The  cracks  located  on  the  axis  of  symmetry  of 
a  fin  were  planar,  grew  in  that  plane  and  exhibited 
Mode  I  fields  all  around  the  flaw  border.  As  such, 
they  were  Class  I  (Cotterell,  1965)  cracks  throughout 
their  growth  and  grew  readily  under  internal  pres¬ 
sure  revealing  no  unusual  effects.  On  the  other  hand, 
the  cracks  which  were  initiated  at  the  confluence  of 
the  R1.3  and  Rll  values  (Fig.  1)  were  initially  non- 
planar  and  proceeded  to  turn  under  internal  pres¬ 
sure,  exhibiting  Class  II  crack  behavior  until  .find¬ 
ing  their  Class  I  direction  for  further  growth.  Crack 
growth  was  retarded  by  this  turning  effect  which  in¬ 
volved  the  elimination  of  shear  modes  present  dur¬ 
ing  turning.  This  study  focuses  on  two  types  of  such 
cracks  called  off-axis  cracks. 

Two  kinds  of  off-axis  cracks  were  used,  those  en¬ 
tering  the  model  normal  to  the  fin  surface  and  those 
entering  the  model  parallel  to  the  fin  axis.  The  first 
type  was  made  by  drilling  a  small  hole  opposite  the 
fin  to  be  cracked,  (Fig.  2)  and  inserting  a  shaft  with 
two  blades  such  that  the  blade  producing  the  crack 
was  normal  to  the  fin  surface.  Typical  cracks  formed 
in  this  way  with  the  blade  held  normal  to  the  surface 
are  shown  in  Fig.  3. 

In  Model  4,  all  of  the  crack  growth  was  under 
mixed  mode  (Modes  I  &  II)  loading  and  both  modes 
existed  everywhere  when  growth  was  stopped  except 
near  the  fin  surface  where  only  Mode  I  existed.  The 
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crack  is  thus  still  a  Class  II  crack.  In  Model  8-i,  the 
crack  has  grown  out  of  both  Modes  II  and  III,  be¬ 
coming  parallel  to  the  fin  axis,  and  exhibited  only  a 
Mode  I  SIF  along  the  crack  border.  The  river  marks 
indicated  that  Mode  III  was  present  earlier.  They 
are  exacerbated  by  slight  tilting  of  the  blade  with 
respect  to  the  fin  surface.  Retardation  of  growth 
along  the  crack  front  by  Mode  III  is  shown  in  Model 
8  in  the  region  of  the  river  marks. 

Fig.  4  shows  a  typical  starter  crack  and  its 
growth  when  the  blade  is  held  parallel  to  the  fin 
axis.  A  small  damaged  region  is  shown  at  entry  into 
the  model  but  the  crack  quickly  orients  itself  parallel 
to  the  fin  axis  under  internal  pressure  and  exhibited 
pure  Mode  I  all  around  the  crack  front  as  a  Class 
I  crack.  A  few  river  marks  indicate  the  presence  of 
some  Mode  III  which  was  eliminated  during  growth. 

It  appeared  that  large  displacements  normal  to 
the  crack  surface  were  responsible  for  the  river  marks 
so  no  calculations  were  made  where  they  were  near 
the  surface  in  any  of  the  cracks. 

Since  turning  occurs  above  the  critical  tempera¬ 
ture  (Tc)  of  the  material,  the  turning  effect  is  grad¬ 
ual  (Smith,  Constantinescu  and  Liu,  2001),  and  not 
a  kink. 

In  examining  the  results,  it  appears  that  the  pres¬ 
ence  of  the  river  markings  exerted  a  substantial  ef¬ 
fect  upon  the  SIF  calculations  when  the  data  were 
taken  before  significant  crack  growth  has  occurred 
after  turning.  This  was  due  to  several  factors: 

i.  The  projected  crack  shape  deviated  signifi¬ 
cantly  from  a  semi-ellipse  during  and  just  after 
turning  (Fig.  3,  Model  8-i) 

ii.  The  fracture  surface  was  not  smoothly  contin¬ 
uous. 

iii.  Where  the  river  markings  extended  to  the 
crack  front,  one  expects  large  non-linear  Mode 
III  effects  (Fig.  3  Model  8-i). 

vi.  Crack  growth  beyond  turning  was  not  sufficient 
to  eliminate  the  significance  of  items  i  through 
iii  for  shorter  cracks  and  so  substantial  scatter 
resulted  in  the  SIF  values  in  these  cases,  and 
are  not  reported  herein. 


A  two  parameter  method  developed  earlier 
(Smith  and  Kobayashi,  1993)  was  used  to  convert 
near  tip  stress  fringe  patterns  into  SIF  values.  The 
algorithms  for  Modes  I  and  II  are  given  in  Appendix 
B.  No  Mode  III  calculations  were  made  due  to  the 
large  out  of  surface  deformations.  Normalized  values 
of  the  SIF’s  (ie.  F{)  were  computed  as  indicated  in 
Table  I  assuming  that  the  cracks  were  semi-elliptic 
and  planar.  This  was  nearly  true  for  the  cracks  en¬ 
tering  the  body  parallel  to  the  fin  axis.  However,  for 
cracks  entering  the  body  normal  to  the  fin  surface, 
cracks  were  neither  planar  nor  semi- elliptic  as  noted 
above. 

For  the  cases  where  the  cracks  were  inserted  par¬ 
allel  to  the  fin  axis,  all  of  the  cracks  were  Class  I 
and  pure  Mode  I  and  results  are  plotted  in  Fig.  5. 
These  cracks  exhibited  negligible  turning  and  their 
final  depths  included  much  more  growth  (25  to  75% 
of  the  final  crack  depth)  than  those  initiated  normal 
to  the  fin  surface  at  the  same  off-axis  locations.  The 
data  scatter  are  all  within  the  experimental  accuracy 
of  6%. 

SUMMARY 

A  series  of  frozen  stress  photoelastic  experiments 
were  conducted  on  internally  pressurized  generic  six 
finned  motor  grain  models.  Starter  cracks  were  initi¬ 
ated  by  sharp  blades  which  produced  real  cracks  by 
wedging  the  material  open.  Crack  growth  and  SIF 
values  were  measured  along  the  flaw  borders.  Results 
showed  that  off-axis  cracks  produced  shear  modes 
which  retarded  the  crack  growth  relative  to  symmet¬ 
ric  cracks  which  initiated  as  Class  I  cracks  under  pure 
Mode  I.  After  turning,  the  off-axis  cracks,  initially 
Class  II  cracks,  became  Class  I  cracks  and  grew  par¬ 
allel  to  the  fin  axis.  Results  showed  that  from  a/t 
of  0.20  to  0.68,  the  normalized  SIF  decreased  with 
increasing  depth. 
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APPENDIX  A-  Frozen  Stress 
Photoelasticity 

When  a  transparent  model  is  placed  in  a  circu¬ 
larly  polarized  monochromatic  light  field  and  loaded, 
dark  fringes  will  appear  which  are  proportional  to  the 
applied  load.  These  fringes  are  called  stress  fringes 
or  isochromatics,  and  the  magnitude  of  the  maxi¬ 
mum  in-plane  shear  stress  is  a  constant  along  a  given 
fringe. 

Some  transparent  materials  exhibit  mechanical 
diphase  characteristics  above  a  certain  temperature, 
called  the  critical  temperature  (Tc).  The  material, 
while  still  perfectly  elastic  wall  exhibit  a  fringe  sen¬ 
sitivity  of  about  twenty  times  the  value  obtained  at 
room  temperature,  and  its  modulus  of  elasticity  will 
be  reduced  to  about  one  six-hundredth  of  its  room 
temperature  value.  By  raising  the  model  tempera¬ 
ture  above  Tc,  loading,  and  then  cooling  slowly  to 
room  temperature,  the  stress  fringes  associated  with 
Tc  will  be  retained  when  the  material  is  returned  to 
room  temperature.  Since  the  material  is  so  much 


more  sensitive  to  fringe  generation  above  Tc  than 
at  room  temperature,  fringe  recovery  at  room  tem¬ 
perature  upon  unloading  is  negligible.  The  model 
may  then  be  sliced  without  disturbing  the  “frozen 
in”  fringe  pattern  and  analyzed  as  a  two-dimensional 
model  but  containing  the  three-dimensional  effects. 
In  the  use  of  the  method  to  make  measurements  near 
crack  tips,  due  to  the  need  to  reduce  loads  above 
critical  temperature  to  preclude  large  local  deforma¬ 
tions,  and  the  use  of  thin  slices,  few  stress  fringes  are 
available  by  standard  procedures.  To  overcome  this 
obstacle,  a  refined  polariscope  is  employed  to  allow 
the  tandem  use  of  the  Post  and  Tardy  methods  to 
increase  the  number  of  fringes  available  locally. 

In  fringe  photographs,  integral  fringes  are  dark  in 
a  dark  field  and  bright  in  a  bright  field. 

Appendix  B 

Mode  I  Algorithm 

Beginning  with  the  Griffith-Irwin  Equations,  we 
may  write,  for  Mode  I,  for  the  homogeneous  case, 

K 

®ij  —  “  TT  /yW  "h  j  (l  3  —  ft  jZ) 

(27 XT)  2 

(1) 

where: 

c Tij  are  components  of  stress 

Kx  is  SIF 

r,6  are  measured  from  crack  tip  (Fig.  B-l) 
are  non-singular  stress  components 

Then,  along  9  =  7r/2,  after  truncating  <7ij 


^.max  _  FIl  I  ^_o  _  FI  A  P 
n~  (87rr)  i  (8ttt)  i 


(2) 


where: 

r°  =  f(<?ij)  and  is  constant  over  the  data  range 
Kap  —  apparent  SIF 

r™ax  =  maximum  shear  stress  in  nz  plane 


Normalizing  with  respect  to  <r, 


cr(ira)2  5{-k  a)  2  c  \a/ 
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where  (Fig.  B-l)  a  =  crack  length,  and  a  =  remote 
normal  stress 
Kap 


i.e 


a(7ra)i 


,s-  ]fli 


is  linear. 


extrapolate  to  r  =  0  to  obtain  K*.  Knowing,  K*f 
K2/K1  and  8^,  values  of  K\  and  K2  may  be  deter¬ 
mined  since 


From  the  Stress-Optic  Law,  r™ax  =nf/2t  where, 
n  =  stress  fringe  order, 

/  =  material  fringe  value,  and 
t  =  specimen  (or  slice)  thickness 
then  from  Eq.  2 


Kap  =  r™zax(8irr)*  =  = 

where  Kap  (through  a  measure  of  n)  and  r  become 
the  measured  quantities  from  the  stress  fringe  pat¬ 
tern  at  different  points  in  the  pattern. 

In  the  present  study,  instead  of  normalizing 
K  with  respect  to  a  (no,)1/2,  we  have  selected 
PsfV^lQ  35  the  normalizing  factor  where  ^fQ  is 
an  elliptic  integral  of  the  second  kind  approximated 
here,  as  shown  in  Table  I.  An  example  of  the  deter¬ 
mination  of  F\  in  Table  I  from  test  data  is  given  in 
Fig.  B-2. 

Mixed  Mode  Algorithm 

The  mixed  mode  algorithm  was  developed  (see 
Fig.  B-3)  by  requiring  that: 

Wl (St rrTO)1/8^M.Ki 

(4) 

w'hich  leads  to: 


K*  =  [(. KlSine°m  +  2  K2cos9°m)2 

+  (K2sin9°J2}i  (6) 

Knowing  K*  and  Ki  &  K2  can  be  determined 
from  Eqs.  5  and  6.  Details  are  found  in  Smith  and 
Kobayashi  (1993). 


Fig.  B-l  Near  Tip  Notation  for  Mode  I. 
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Model  8,  off-axis  inclined. 

a/c  =  0.59,  a /t  =  0 34,  depth  (d) 

y  =  !.00x  + 1.99  — - - * 

#*-^*^*  * 

'  7JZ 

By  measuring  8 ^  which  is  approximately  in  the  di¬ 
rection  of  the  applied  load,  K2/K1  can  be  deter¬ 
mined. 

Then  writing  the  stress  optic  law  as: 

max  _  fn  _  Kap 
n 2t  (8?rr)  2 


where  K, 

K 


&{TTa)^ 


ap 

VS. 


is  the  mixed  mode  SIF,  one  may  plot 
‘s/rja  as  before,  locate  a  linear  zone  and 


Fig.  B-2:  Determination  of  F\  for  Test  Data. 


Fig.  B-3:  Determination  of  6^  for  Mixed  Mode 
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Table  I  DATA  &  RESULTS 


Loads1 

Crack  Description2 

(dimensions  in  mm) 

Model  4 

P  =  88.97  N 

Off-axis  inclined 

Pmax  -  0.049  MPa 

a  —  8.71  Aa  =  2. 1 8 

Fj  =  1 .90 

Notations: 

ps/=  0.035  MPa 

c  =  11.15  Ac  =  3.02 
a/c  =  0.78  aft-  0.23 

.F;  =  0.48 

1 .  P  =  axial  compressive  load 

P  =  88.97  N 

Model  8-i 

Pmax  =  internal  pressure  to 

Off-axis  inclined 

grow  crack 

Pmax  —  0. 1 03  MPa 
ps/=  0.049  MPa 

a  =  12.50  Aa  =  3.4 
c  =  21.1  Ac  =  10.4 

1.99 

psf=  stress  freezing  pressure 

a/c  =  0.59  a/t  =  0.34 

2.  a  =  crack  depth; 

Model  6 

A  a  =  crack  growth 

P  =  88.97  N 

Off-axis  straight  in 

1.72 

c  =  half  length  of  crack  in 

Pmax  =  0.049  MPa 

a  =  11.60  Aa  =  4.67 

fin  tip  surface 

ps/=  0.035  MPa 

c=  17.00  Ac  =  10.66 

Ac  =  half  crack  growth  in 

o/c  =  0.68  a/t=  0.31 

fin  tip  surface 

Off-axis  straight  in 
a  =  11.23  Aa  =  5.86 

1.86 

3.  F,  =  X,F>lp,F~a 

c=  13.00  Ac  =  6.65 
a/c  =  0.86  a/t  =0.30 

i  =  1,2  at  maximum  depth 

Model  7 

Off-axis  straight  in 

=  approximation  of 

P  =  88.97  N 

a  =  15.60  Aa  =  10.0 

1.58 

elliptic  integral  of  second 

Pmax~  0.103  MPa 

c  =  26.45  Ac  =17.57 

kind 

ps/=  0.049  MPa 

a/c  =  0.59  a/t  =  0.42 

(  V'65 

Off-axis  straight  in 

£=1  +  1.464  -  -<1 

a  =  13.90  Aa  =  4.05 

1.87 

W  c 

c=  18.65  Ac  =  10.17 
a/c  =  0.74  a/t  =0.37 

All  flaws  were  characterized  as 

Model  8-s 

semi-elliptic  flaws  of  depth  a 

P  =  88.97  N 

Off-axis  straight  in 

and  length  2c. 

Pmax-  0.103  MPa 

a  =  7.90  Aa  =  2.8 

1.93 

However,  off-axis  cracks  were 

psf-  0.049  MPa 

c=  13.35  Ac  =  7.75 

neither  perfectly  semi-elliptic 

a/c  =  0.59  a/t  =  0.21 

nor  planar. 

Model  9 

P  =  88.97  N 

Off-axis  straight  in 

Pmax  =  0.1 03  MPa 

a  =  25.10  Aa=  18.7 

1.50 

psf-  0.049  MPa 

c  =  39.4  Ac  =  33.8 
a/c  =  0.64  a/t  =  0.68 
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all  dimensions  are  in  mm 


Fig.  1 :  Dimensions  and  Crack  Locations  in  a  Cross-Section  of  the  Model 


Fig.  2:  Setup  for  Producing  Off-Axis  Starter  Crack 
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(  s  Section  S-S 


Starter  crack  Tip  of  starter  crack 


FS  -  fin  surface  magnification  factor:  2.25 

C  -  crack  front 

D  -  camera  view  of  the  photograph 

Model  7  Off-Axis  Straight  Crack  Showing  Starter  Crack  and  Final  Crack  Front 
Fig.  4:  Typical  Off-Axis  Straight-in  Crack  (Blade  Oriented  Parallel  to  Fin  Axis) 


maximum  depth  normalized  experimental  SIF 


Fig.  5:  Effect  of  Part-Through  Crack  Depth  on  Normalized  SIF 
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